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ABSTRACT 

The relative population of the fine structure sublevels of an atom's ground state is affected by 
radiative transitions induced by an anisotropic radiation flux. This causes the alignment of atomic 
angular momentum. In terms of observational consequences for the interstellar and intergalactic 
\Q , medium, this results in the polarization of the absorption lines. In the paper we consider the conditions 

necessary for this effect and provide calculations of polarization from a few astrophysically important 
atoms and ions with multiple upper and lower levels for an arbitrary orientation of magnetic fields to 
s ! . the a) source of optical pumping, b) direction of observation, c) absorbed source. We also consider an 

astrophysically important "degenerate" case when the source of optical pumping coincides with the 
source of the absorbed radiation. We present analytical expressions that relate the degree of linear 
polarization and the intensity of absorption to the 3D orientation of the magnetic field with respect 
to the pumping source, the source of the absorbed radiation, and the direction of observations. We 
discuss how all these parameters can be determined via simultaneous observations of several absorption 
lines and suggest graphical means that are helpful in practical data interpretation. We prove that 
studies of absorption line polarization provide a unique tool to study 3D magnetic field topology in 
various astrophysical conditions. 

qq ' Subject headings: ISM: atomic processes — magnetic fields — polarization 
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1. INTRODUCTION 



Magnetic fields play extremely important roles in many astrophysical circumstances, e.g., the interstellar medium, 
intergalactic medium and quasars, etc. Unfortunately, there are only a few techniques for magnetic field studies that 
i-Q ' are available (see below). Each technique is sensitive to magnetic fields in particular environments. Therefore even 
Q-f the directions of magnetic field obtained for the same region of sky with different techniques differ substantially. The 
q . simultaneous use of different techniques provides a possibility of magnetic field tomography. 

Polarimetry of aligned dust provides a way of studying magnetic field direction in the diffuse interstellar medium, 
\ molecular clouds, circumstellar and interplanetary medium (see review by Lazarian 2003). Substantial progress in 
d • the understanding of grain alignment has been achieved in the last decade making the technique more reliable. The 
technique is, however, challenging to apply to low column densities, as the polarization signal becomes too weak. It 
. £h '. nas other limitations when dealing with high densities (Cho & Lazarian 2005). 

Polarimetry of some molecular lines using the Goldreich-Kylafis (1982) effect has recently been shown to be a good 
tool for magnetic field studies in molecular clouds (see Girart, Crutcher & Rao 1999). However, the magnetic field 
5^ \ direction obtained has an uncertainty of 90 degrees, which may be confusing. The technique is most promising for 
dense CO clouds. 

Zeeman splitting (see Crutcher 2004) provides a good way to get magnetic field strength, but the measurements are 
very time consuming and only the strongest magnetic fields are detectable this way (see Heiles & Troland 2004). 

Synchrotron emission/polarization as well as Faraday rotation provide an important means to study magnetic fields 
either in distinct regions with strong magnetic fields or over wide expanse of the magnetized diffuse media (Haverkorn 
2005). 

Here we discuss yet another promising technique to study magnetic fields that employs optical and UV polarimetry. 
This technique can be used for interstellar 3 , and intergalactic studies as well as for studies of magnetic fields in QSOs 
and other astrophysical objects. As we discuss below, the technique is based on the ability of atoms and ions to be 
aligned by external radiation in their ground state and be realigned through precession in magnetic field. Henceforth, 
we shall not distinguish atoms and ions and use word "atoms" dealing with both species. This is justifiable, as the 
alignment of ions and atoms does not differ. 

The physics of atomic alignment is rather straightforward. It has been known that atoms can be aligned through 
interactions with the anisotropic flux of resonance emission (see review Happcr 1972 and references therein). Alignment 
is understood here in terms of orientation of the angular momentum vector J, if we use the language of classical 
mechanics. In quantum terms this means a difference in the population of sublevels corresponding to projections of 
angular momentum to the quantization axis. Later we shall refer to J alignment in the ground state of atoms as 
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the "atomic alignment". Note that we deal with the ground state as it is long lived and therefore sensitive to weak 
magnetic field. 

It is worth mentioning that atomic alignment was studied in laboratory in relation with early day maser research 
(see Hawkins 1955). This effect was noticed and employed in the interstellar case by Varshalovich (1968) in the case 
of hypcrfine structure of the ground state. Varshalovich (1971) suggested that atomic alignment enables us to detect 
the direction of magnetic fields in the interstellar medium, but did not provide a thorough quantitative study. 

A more rigorous study of using atomic alignment to diagnose weak magnetic field in diffuse media was conducted 
in Landolfi & Landi Degl'Innocenti (1986). However, in that case, an idealized two-level atom was considered. In 
addition, polarization of emission from this atom was discussed for a very restricted geometry of observations, namely, 
the magnetic field is along the line of sight and both of these directions are perpendicular to the incident light. This 
made it rather difficult to use this study as a tool for practical mapping of magnetic fields in various astrophysical 
environments. 

In this paper, we shall consider the alignment of atoms with fine structure of the ground level. Realistic atomic 
species will be studied, taking into account their multi-level structure. We shall show that atomic alignment can 
influence the line polarization and intensity of absorption lines. The polarization dependence on the 3D geometry of 
both radiation and magnetic field will be provided. A study of the alignment of species within hyperfine structure of 
the ground state (e.g. Na, HI) and the polarization of emission lines will be provided in our companion paper. 

It should be pointed out that the atomic alignment we deal with in this paper differs from the Hanle effect that 
solar researcher have studied. Hanle effect is depolarization and rotation of the polarization vector of the resonance 
scattered lines in the presence of a magnetic field, which happens when the magnetic splitting becomes comparable to 
the decay rate of the excited state of an atom. The research into emission line polarimetry resulted in important change 
of the views on solar chromosphere (see Landi Degl'Innocenti 1983a, 1984, 1998, Stenflo & Keller 1997, Trujillo Bucno 
& Landi Degl'Innocenti 1997, Trujillo Bueno 1999, Trujillo Bueno et al. 2002, Manso Sainz & Trujillo Bueno 2003). 
However, these studies correspond to a setting different from the one we consider here. In this paper we concentrate 
on the weak field regime, in which it is the atoms at ground level that are repopulated due to magnetic precession, 
while the Hanle effect is negligible for the upper state. This is the case, for instance, of the interstellar medium. The 
polarization of absorption lines that we study here is thus more informative. In many cases, we are in optically thin 
regime, so we do not need to be concerned about radiative transfer. 

In what follows we present the discussion of basic atomic alignment physics and optical pumping in §2. Detailed 
treatment of statistical equilibrium of fine structure atoms in a magnetic field is given in §3. Following this we then 
calculate, for a number of elements, their alignment and the induced polarizations of their emission and absorption 
lines in §4. To help our reader understand the logic of the calculations we provide a quantitative level summary of the 
procedures in §5. The discussion and the summary are provided in, respectively, §6 and §7. The Appendix contains 
the information necessary for understanding our formalism as well as auxiliary calculations. 

2. CONDITIONS FOR ATOMIC ALIGNMENT 

The basic idea of the atomic alignment is quite simple. The alignment is caused by the anisotropic deposition of 
angular momentum from photons 4 . In typical astrophysical situations the radiation flux is anisotropic (see Fig[T]). As 
the photon spin is along the direction of its propagation, we expect that atoms scattering the radiation from a light 
beam to be aligned. Such an alignment happens in terms of the projections of angular momentum to the direction of 
the incoming light. For atoms to be aligned, their ground state should have non-zero angular momentum. Therefore 
fine (or hyperfine) structure is necessary to enable various projection of atomic angular momentum to exist in their 
ground state. 

2.1. Basics of atomic alignment 

Let us discuss a toy model that provides an intuitive insight into the physics of atomic alignment. Consider an atom 
with its ground state corresponding to the total angular momentum 1=1 and the upper state corresponding to the 
angular momentum 1 = (Varshalovich 1971). If the projection of the angular momentum to the direction of the 
incident resonance photon beam is M, the upper state M can have values —1, 0, and 1, while for the upper state 
M=0 (see Figl^e/t). The unpolarized beam contains an equal number of left and right circularly polarized photons 
whose projections on the beam direction are 1 and -1. Thus absorption of the photons will induce transitions from the 
M = — 1 and M = 1 sublevels. However, the decay from the upper state populates all the three sublevels on ground 
state. As the result the atoms accumulate in the M = ground sublevel from which no excitations are possible. 
Accordingly, the optical properties of the media (e.g. absorption) would change. 

The above toy model can also exemplify the role of collisions and magnetic field. Without collisions one may expect 
that all atoms reside eventually at the sublevel of M = 0. Collisions, however, redistribute atoms to different sublevels. 
Nevertheless, as disalignment of the ground state requires spin flips, it is less efficient than one might naively imagine 
(Hawkins 1955). The reduced sensitivity of aligned atoms to disorienting collisions makes the effect important for 
various astrophysical environments. 

A magnetic field would mix up different M states. However, it is clear that the randomization in this situation is 
not complete and the residual alignment reflects the magnetic field direction (see Fig[T]). Magnetic mixing happens if 

4 We do not consider interactions with other particles (see Fineschi & Landi Degl'Innocenti 1992). 



3 




Fig. 1. — Typical astrophysical environment where atomic alignment can happen. A pumping source deposits angular momentum to 
atoms in the direction of radiation and causes differential occupations on their ground states. In a magnetized medium where the Larmor 
precession rate vl is larger than the photon arrival rate , however, atoms are realigned with respect to magnetic field. Atomic alignment 
is then determined by 9 r , the angle between the magnetic field and the pumping source. Observed polarization depends on both 8 r and 
9, the angle between the magnetic field and the line of sight. In general, there are two situations: left, the alignment is produced by a 
pumping source while we observe another weak background source whose light passes through the aligned medium; right, the background 
source coincides with the pumping source, in this case, 9 r = 9. 




Fig. 2. — Left: a toy model to illustrate how atoms are aligned by anisotropic light. Atoms accumulate in the ground sublevel M = as 
radiation removes atoms from the ground states M = 1 and M = — 1; Right: Radiation geometry and the polarization vectors in a given 
coordinate system. O is the direction of radiation, \ is the positional angle of a linear polarization. 

the angular momentum precession rate is higher than the rate of the excitation from the ground state, which is true 
for many astrophysical conditions, e.g., interplanetary medium, ISM, intergalactic medium, etc. 

Since conservation and transfer of angular momentum is essential for the problem we deal with, it is necessary to 
quantize the radiation field as the atomic states are quantized. The classical theory can give a qualitative interpretation 
which shall be utilized in this paper to provide an intuitive picture. However, quantitative results have to be obtained 
by the quantum- mechanical treatment, namely, density matrix description for both atoms and the radiation. The key 
issue here is that the focus is on angular momentum instead of energy. As we are interested in exchange of angular 
momentum between photons and atoms, we have to quantize the radiation field. 
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Fig. 3. — Different regimes divided according to the strength of magnetic field and the Doppler line width. Atomic realignment is 
applicable to weak field (< 1G) in diffuse medium. Level interferences are negligible unless the medium is substantially turbulent (Sv > 
lOOkm/s) and the corresponding Doppler line width becomes comparable to the fine level splitting u&j. For strong magnetic field, Zeeman 
effect dominates. When magnetic splitting becomes comparable to the Doppler width, a and 7r components can still distinguish themselves 
through polarization, this is the magnetograph regime; Hanle effect is dominant if Larmor period is comparable to the lifetime of excited 
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similarly, for ground Hanle effect, it requires Larmor splitting to be of the order of photon pumping rate; for weak 



magnetic field (< 1G) in diffuse medium, however, atomic alignment is the main effect provided that i>l > 
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TABLE 1 

Selected alignable atomic s pecies and corresponding transitions. No te only lines above 912A are listed. Data are taken 
from the Atomic Line List |http://www.pa.uky.edu/~peter/atomic/| and the NIST Atomic Spectra Database. The last 

column gives the maximum polarizations and its corresponding transitions. For those species with multiple lower levels, 
the polarizations are calculated for shell star (t e ff = 15, 000k) in the strong pumping regime; in the weak pumping 
regime, the maximum polarizations are 19% for oi transition (2 — > 2), and 9% for si transition (2 — > 2). 



All in all, in order to be aligned, first, atoms should have enough degrees of freedom: namely, the quantum angular 
momentum number must be > 1. Second, the incident flux must be anisotropic. Moreover, the collisional rate 
should not be too high. While the latter requires special laboratory conditions, it is applicable to many astrophysical 
environments such as the outer layers of stellar atmospheres, the interplanetary, interstellar, and intergalaetic medium. 

2.2. Relevant time-scales 

Various species with fine structure can be aligned. A number of selected transitions that can be used for studies of 
magnetic fields are listed in Table [TJ Why and how are these lines chosen? First, we gathered all of the prominent 
interstellar and intergalaetic lines (Morton 1975, Savage et al. 2005), from which we take only alignable lines, namely, 
lines with ground angular momentum number J g > 1. Then we rule out those species with nuclear spin since hyperfine 
lines will be dealt with in our companion paper. The number of prospective transitions increases considerably if we 
add QSO lines. In fact, many of the species listed in the Table 1 in Verner, Barthel, & Tytler (1994). are alignable and 
observable from the ground because of the cosmological redshifts. In this paper we do not deal with such transitions. 

In terms of practical magnetic field studies, the variety of available species is important in many aspects. One of 
them is a possibility of getting additional information about environments. Let us illustrate this by considering the 
various rates (see Table [5]) involved. Those are 1) the rate of the Larmor precession, z/£, 2) the rate of the optical 
pumping, Tjj 1 , 3) the rate of collisional randomization, t^ 1 , 4) the rate of the transition within ground state, 1 . 
In many cases vl > t^ 1 > t~ 1 ,t^I . Other relations are possible, however. If r^ 1 > t^ 1 , the transitions within the 
sublevels of ground state need to be taken into account and relative distribution among them will be modified (see 
i )4.2.2p . Since emission is spherically symmetric, the angular momentum in the atomic system is preserved and thus 
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TABLE 2 

Examples of absorptions from metastable states. Calculation for Cr II shows that the absorption from metastable level 
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TABLE 3 

Relevant rates for atomic alignment. A m is the magnetic dipole emission rate for transitions among J levels of the 

GROUND STATE OF AN ATOM. f k j IS THE INELASTIC COLLISIONAL TRANSITION RATES WITHIN GROUND STATE DUE TO COLLISIONS WITH 
ELECTRONS OR HYDROGENS, AND f ap IS THE SPIN FLIP RATE DUE TO VAN DER WAALS COLLISIONS (SEE i|4.2.3t . In THE LAST ROW, EXAMPLE 
VALUES FOR C II ARE GIVEN. T^ 1 IS CALCULATED FOR AN O TYPE STAR, WHERE R* IS THE RADIUS OF THE STAR RADIUS AND R IS THE 

DISTANCE TO THE STAR. 

alignment persists in this case. In the case vl < r# , the magnetic field does not affect the atomic occupations and 
atoms are aligned with respect to the direction of radiation (see FigHJ. From the expressions in Table [31 we see, for 
instance, that magnetic field can realign CII only at a distance r > 7.7Au from an O star if the magnetic field strength 
~5/iG. 

If the Larmor precession rate vl is comparable to any of the other rates, the atomic line polarization becomes sensitive 
to the strength of the magnetic field. In these situations, it is possible to get information about the magnitude of 
magnetic field. In this paper we consider the regime where vl is much smaller that the decay rate of the excited 
state, which means that we disregard the Hanle effect. Fig[3] illustrates the regime of magnetic field strength where 
atomic realignment applies. Atoms are aligned by the anisotropic radiation at a rate of . Magnetic precession will 
realign the atoms in their ground state if the Larmor precession rate i>l > r^ 1 . In contrast, if the magnetic field gets 
stronger so that Larmor frequency becomes comparable to the line- width of the upper level, the upper level occupation, 
especially coherence is modified directly by magnetic field, this is the domain of Hanle effect, which has been extensively 
discussed for studies of solar magnetic field (see Landi Degl'Innocenti 2004 and references therein). When the magnetic 
splitting becomes comparable to the Doppler line width v>jj, polarization appears, this is the "magnetograph regime" 
(Landi Degl'Innocenti 1983b). For magnetic splitting vl ^> vb-, the energy separation is enough to be resolved, and 
the magnetic field can be deduced directly from line splitting in this case. If the medium is strongly turbulent with 
Sv ~ lOOkm/s (so that the Doppler line width is comparable to the level separations i>d ~ ^aj), interferences occur 
among these levels and should be taken into account. 

Long-lived alignable metastable states that are present for some atomic species between upper and lower states may 
act as proxies of ground states. Absorptions from these metastable levels can also be used as diagnostics for magnetic 
field therefore. We present a few selected species with such states in Table [2 but further in the present paper deal only 
with one representative of this class, namely, Cr II. The life time of the metastable level is an important timescale. 
However, in the formalism of the present paper one can account for it via modifying the effective tt to account for the 
decay of the metastable states. 

2.3. Ground state alignment and absorptions 

The geometry of the radiation system is illustrated by FigHJe/t. The origin of this frame is defined as the location 
of the atomic cloud. The line of sight defines z axis, and together with direction of radiation, they specify x-z plane. 
The x-y plane is thus the plane of sky. In this frame, the incident radiation is coming from (6q, 0), and the magnetic 
field is in the direction (9b,4>b)- 

The magnetic field is chosen as the quantization axis (z") for the atoms. Alignment shall be treated in the frame 
x^y^z" (see ApplElhow these two frames are related). In this "theoretical" frame, the line of sight is in (9, tt) direction 
(i.e., the x"-z" plane is defined by the magnetic field and the line of sight, see Fig\4}right) , and the radiation source is 
directed along (6 r ,<fr r ). 

Absorption is directly associated with the alignment of the ground state and the direction of observation (Fig[T]). For 
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Fig. 4. — Left: geometry of the observational frame. In this frame, the line of sight is z axis, together with the incident light, they 
specify the x-z plane. Magnetic field is in (0b,</>b) direction; Right: transformation to the "theoretical frame" where magnetic field defines 
z" axis. This can be done by two successive rotations specified by Euler angles (4>B,9r) (see App[E]for details). The first rotation is from 
xyz coordinate system to x'y'z' coordinate system by an angle 4>B about the z-axis, the second is from x'y'z' coordinate system to x"y"z" 
coordinate system by an angle about the y'-axis. Atomic alignment and transitions are treated in the "theoretical" frame where the line 
of sight is in (9,ir) direction and the incident radiation is in (6 r ,ipr) direction. 

absorption, we shall demonstrate later that the Stokes parameter (see App|D]) U = if the background is unpolarized. 
This means that the linear polarization is either parallel or perpendicular to the magnetic field. This happens because 
of fast precession of atoms in the ground state. Direction of polarization thus traces 4>b with a 90° degeneracy. For 
absorptions in the optically thin case 5 , the Stokes parameters Q = —r]idIo,I = Iq(1 — r/od), where rji are absorption 
coefficients and given by Eq lf 41 d is the thickness of the medium. Iq is the intensity of background source. Thus the 
degree of polarization P per unit optical depth t = r/ad is given by 

P = _Q_ „ zJh = 1.5<7g(J,)sin 2 0w 2 JlJu 

t Irj d~ t] s/2 + ul{Ji){l-l.hsm 2 e)w 2 j iJ J 

<7q depends on 9 r and is the normalized dipole component of the ground state density matrix, which is the main 
parameter we calculate in this paper, w^j is defined by Eq.(|13p and the corresponding numerical values are given in 
Table [2] For a generic case, where the background source (e.g. QSO) is polarized and optical depth is finite 6 , we can 
get in the first order approximation 

/ = (/n + Qo)e~ T{1+rllM + (I - Qo)e- r(1 -" l/?)o) , 
Q = (Io + g )e- T(1+ " l/ " o) - (Jo - Qo)e- T{1 -^ \ 

U = U Q e~ T 7 V = Voe~ T , (2) 

in which Iq,Q q ,Uq, Vq are the Stokes parameters of background source. Optical depth is also modified by the alignment. 
The line ratio of a doublet varies as 

U = i?./ iJl [V2 + f 7 2 (J ; )(l - 1.5sin 2 0)w 2 JlJ2 ] 
t 2 B Jl j 2 [V2 + al(J l ){l-l.hs^e)w\ j y 

where Bj t j u represents Einstein coefficient. The next section will be devoted to explaining how the atomic alignment 
(<7q) can be calculated using QED theory. Readers who are not interested in the technical details, however, can skip 
that section and go directly to the results of calculations for the species in Table 1 that we present in SjU The value 
of fj3] is that it provides the approach for calculating the alignment of atoms other than those discussed in the present 
paper. 

3. ALIGNMENT THEORY 

In what follows, we consider a realistic atomic system, which can have multiple upper levels J U1 M u and lower levels 
J/,M;. When such an atomic system interacts with resonant radiation, there will be a photoexcitation followed by 
spontaneous emission. The atomic occupation is thus determined by the balance of the three processes: absorption, 
emission and magnetic precession among the sublevels of the state. The radiative transition rates can be constructed 

5 In solar case, one has to deal with radiative transfer. For diffuse medium, however, this assumption is reasonable. 

6 We neglect emission here. 
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from the transition probabilities given by Eq. (|Bl[) . Then additional transformation should be made to express the 
density matrices in irreducible form (see App[A}. The actual derivation of the equations that govern the statistical 
equilibrium of the atomic occupations is rather complicated and can be found in Landi Dcgllnnocenti & Landol (2004). 
Below are the statistical equilibrium equations for the upper and ground state (Degllnnocenti & Landol 2004): 



P k q {ju)- 



-i2TTV L g u qp k {J u ) = 



(3[k,k',K 




A{J u ^Ji)p k q (J u )+ [ J iWi 

JiKQk'q' 



Ju){-lf +q ' 



kk'K 
qq'Q 



jQP- q '(Ji) 



(4) 



p%(J l )+i2nv L g l qp k (J l ) = J^(-l) J ' +J - +k+1 [Ju] {jJll} MJu - Ji)p k (J u ) 
-[Jl] J2 B{J l ^J u )^l)' h - J - +k+k '+'i' +1 {-i[k,k',K]) 1 ! 2 

J^KQk'q' 

1 I K^ { k k'K) fkk' K\ JK k > / t \ 



Jl Jl Ji, 



Jl Jl Jl 



qq'Q 



(5) 



where p q and Jq are the irreducible form of density matrix of atoms and radiation, respectively. The evolution of 
upper state (Pg(J u )) is represented by Eq.Q and the ground state (p k {Ji)) is described by Eq.([5]). The second terms 
on the left side of Eq. (|4|5p represent mixing by magnetic field, where g u and gi are the Lande factors for the upper 
and ground level. For the upper level, this term can be ignored as we consider a regime where the magnetic field is 
much weaker than Hanle field 7 (vl <C A/g u ). The two terms on the right side of Eq.(|U [5]) are due to spontaneous 
emissions and the excitations from ground level. Transitions to all upper states are taken into account by summing 
over J u in Eq.©. Vice versa, for an upper level, transitions to all ground sublevels (Ji) are summed up in Eq.(j4]). 
Let us remind our reader, this multi-level treatment is well justified as level crossing interferences do not exist in the 
weak field regime unless the medium gets strongly turbulent so that Doppler width becomes comparable to the energy 
separation between levels (Sv > lOOkm/s, see FigEJ). The matrix with big "{ }" represents the 6-j or 9-j symbol, 
depending on the size of the matrix. When more than two angular momentum vectors are coupled, there is more than 
one way to add them up and form the same resultant. The 6-j (or 9-j) symbol appears in this case as a recoupling 
coefficient describing transformations between different coupling schemes (see ApplC]). Throughout this paper, we 
define [j] = 2j + 1, which means [Ji] = 2Ji + 1 and [J u ] = 2J U + 1. The Einstein spontaneous emission rate is defined 
as the total probability of an atom from a specific state J', M' to all states of J: A = 
for stimulated transitions B are related to A by: 



M ' 



The Einstein coefficients 



[Ju}A(J u Ji) = 2hv 3 /c 2 [J l ]B(J l J u ) = 2hv z /c 2 [J u ]B{J u -+ Ji), 



(6) 



where [J] is the statistical weight of the initial level J. Stimulated emission is neglected since A/g u ^S> vl > t^ 1 ~ 
B(J U — > Ji)Jq- Note for spontaneous emission and magnetic mixing, k,q are conserved because of the symmetric 
nature of these processes. Indeed, magnetic mixing and emission are axially symmetric processes. 
The excitation depends on 



Jq = / dv^i(v - v ) 



4ir ^ 

i=0 



(7) 



which is the radiation tensor of the incoming light averaged over the whole solid-angle and the line profile £(y — vq). 
Si = [I, Q, U, V] represent the Stokes parameters (App|D]). In many cases, the radiation source is far enough that it 
can be treated as a point source. We consider here unpolarized incident light, thus Q = U = V = 0. The irreducible 
unit tensors for Stokes parameters /, Q, U fApp lDj) are: 



J$(i,n)= o , j 2 (i,n) = V2 



(1 - 1.5 sin 2 9)/2 
-3/4 sin 2 (9 




Jl 2 (i,n) = Vs 



sin 2 0/4 



-(1 + cos 2 0)/4 
=F« cos 0/2 



Tsin20/4\ 
J±i(*,fi) = %/3 ( T sin 20/4 
i sin 0/2 / 



,±id 



(8) 



They are determined by its direction f2 and the reference chosen to measure the polarization (¥\g^right). As addressed 
earlier, the line of sight is in Q — (0, tt) in the "theoretical frame" . For the incoming radiation from (0 r , <f> r ), according 
to Eq.®, the nonzero elements of the radiation tensor are: 



7 For the Hanle effect to be dominant, magnetic splitting ought to be comparable to the energy width of the excited level; see Fig|3] 
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J' 


3/2 5/2 7/2 
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0.5292 -0.6047 0.1890 
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0.4629 -0.6172 0.2160 


J 


4 


9/2 


5 


J' 


3 4 5 


7/2 9/2 11/2 


4 5 6 




0.4432 -0.6205 0.2256 


0.4282 -0.6228 0.2335 


0.4163 -0.6245 0.2402 



TABLE 4 

Numerical values of w 2 , ,,. J is the J value of the initial level and J' is that of the final level. 



70 r 72 W a to o - 2 n \ j 72 VSWg 



sin 2 9 r he ±2i4 



f2 _ VZWg 



sw.29 r he ±i4 



(9) 



where W is the dilution factor of the radiation field, which can be divided into anisotropic part W a and isotropic part 
Wi (Bommier & Sahal-Brechot 1978). In the case of a point source, Wi = 0. If Wi 7^ 0, the degree of alignment and 
polarization will be reduced. The solid-angle averaged intensity for a black-body radiation source is 

h = w^- — L r . (10) 



e hv/k B T _ l ■ 

By setting the first terms on the left side of Eq. (|4l5|) to zeros, we can obtain the following linear equations for the 
steady state occupations in the ground states of atoms: 

AlJ " h) 7T Y} J iWi - J «) E p-»> ( J /) E(- 1 ) fc+fe ' +9 ' 



2-nip k q {Ji)qgiv L - ^ 



A{J U -> J! 



k'q' 



KQ 



(3[k,k',K])V 2 



(-1) J '- J "6., 



kk'K 
qq'Q 

1 1 K 



IK 



(_1)J.+J«[J, 



l JuJuL ) I k k W 



JiJ n Ji Ji j u i \ Ji Ji Ji 



k k' K 







(11) 



The multiplet effect is counted for by summing over J u . The alignment of atoms induces polarization for absorption 
lines even though the background radiation is unpolarized. The corresponding absorption coefficients rji are (Landi 
DeglTnnocenti 1984) 

Tfcfofi) = ^Bn^JJ^-^B-llMj^W)^^"). ( 12 ) 

n KQ 



where n(Ji) = n\/[Ji]Po(Ji) is the total atomic population on level J/, (Tq = Pq / Pq, and 



4j„ 



1 1 K 

Jl Jl Ju 



1 1 

Jl Jl Ju 



(13) 



The values of Wjj, for different pairs of J, J' are listed in Table [4j 
Since there is no magnetic coherence because of fast magnetic mixing on ground level, cr^ = 0, 



(14) 



K 



From the above equation, we see that r/2 = as J${2,tt) = (EqJSJ) for our choice of the reference direction of 
polarization. In this case x is either 0° or 90°, which means the polarization is either parallel or perpendicular to the 
projection of magnetic field in the plane of sky (Fig^-ight) . 

Comparing with absorption, emission is negligible as emitted intensity I e m oc p u oc Jg (the solid-angle averaged 
intensity) which is reduced by a dilution factor W 1 fEa lTO]) , 

Thereby for optically thin case, the Stokes parameters of an absorbed unpolarized radiation are, 



I = (1 - r) d)Io, Q = -Vidlo 



(15) 



If we drop the assumption of unpolarized background source, the more general expression for finite optical depth would 
be 

I=(Io + Qo)e- ( " 0+r ' l)d + (Jo - Q )e-^-^ d , Q = (I + Q )e- ( " 0+r ' l)d - (J - Qo)?^"^, U = U^-\ V = V e~ T 

(16) 
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Fig. 5. — Left: The normalized density tensor components <jg = p\l p% of the ground state of S II, dash-dot line is the result if we only 
consider the transition (45°y 2 — 4P 3 / 2 ) itself; solid line is the actual result with the multiplet effect taken into account, the transitions 

between the ground state 4S'°y 2 and other upper levels 4P 3 y 2 ,5/2 reduce the degree of alignment; Right: The normalized dipole density 

components Oq = Pq/Pq of the ground state aGS ('x' line) and metastable state 16.03/2, 5/2, ...9/2 (other lines) of Cr II. 

where Iq, Qq, Uq are the Stokes parameters of the background radiation, which can be from a weak background source 
or the pumping source itself, d refers to the thickness of medium. 

For unpolarized pumping light, the radiation has only components with K = 0, 2. As the result, the density tensor 
of the atom has only even k components, while odd k components will vanish according to the selection rules of "3j", 
"6j", "9j" symbols (see App. [UJ. Physically, this means atoms are aligned but not oriented. 

4. ALIGNMENT OF VARIOUS ATOMIC SPECIES 

4.1. Single ground level species: S II, & Cr I 

The S II and Cr I absorption lines are observed in interstellar medium (Morton 1975). We shall discuss S II first. 
S II has a ground level ^S^ 2 and upper levels 4Pi/2,3/2,5/2) which means J/ = 3/2, J u = 1/2,3/2,5/2. According to 

the triangle rule (AppJC]) of the 3j symbol in Eq.jATJ, the irreducible density tensor of the ground state p\{Ji) has 
components with k = 3/2 — 3/2, 3/2 — 1/2, 3/2 + 3/2 = 0, 1, 2, 3. For unpolarized pumping light, we only need to 
consider even components k = 0, 2 as explained earlier. 

The fact that vl 3> t^ 1 indicates that the coefficients of the diagonal terms with q ^ in Eq. fTT]) are much larger 
than the other coefficients. This implies all Zeeman coherence components (q ^ 0) disappear. It is expected that 
there is no Zeeman coherence in the magnetic quantization frame as magnetic mixing is fast. Owing to the selection 
rule of "3j" symbols (see App[C|), only Jq'_q components appear and they are determined by the polar angle 6 r of the 

radiation (Eq|5J FigU]). As a result, Pq' 4 are real quantities and they are independent of the azimuthal angle 4> r (see 
Eq[5]). Physically this results from fast procession around magnetic field. Their dependence on polar angle is shown in 
Fig[5j To calculate the alignment for the multiplet, all transitions should be counted according to their probabilities 
even if one is interested in only one particular line. This is because they all affect the ground populations and therefore 
the degree of alignment and polarization. This means that in practical calculations, a summation should be taken 
over all the three upper levels J u = 1/2, 3/2, 5/2 in Eq.fTT]). Inserting the values of J u , Ji, k, k' , K and Jq (EqJ5]) into 



Eq.fTJJ), we get three linear equations for p ' (J;) 



By solving them, we obtain 

1.0481 -3. 1442 cos 2 r 
0.264 cos 2 r - 13.161 



(17) 



where the floating point numbers are approximate. The same is true with other results in the paper; we take > 4 
significant digits. To demonstrate the effect of multiplet, we also calculate the alignment resulting from only one 
transition 3/2 — ► 1/2. The comparison is given in FiglSl Apparently, multiple transitions reduce the degree of 
alignment. 

Substituting Eq. (jTTJ) into Eq. (jT|) , we get the polarization of S II absorption line, 



(1.5722 - 4.7163 cos 2 9 r ) sin 2 6wj j 



(1.4822 - 4.4466 cos 2 8 r ) - (1.0481 - 3. 1442 cos 2 r )(l - 1.5sin 2 



(18) 



FigEl shows the dependence of the degree of polarization on 8 r and 9. FigJ5]!e/£ is the corresponding contour plot. We 
see that the polarization in absorption is either parallel or perpendicular to the plane of sky. The polarization changes 
sign at the Van Vlcck angle 54.7° (Van Vleck 1925, House 1974). This happens because of coupling of the two oscillators 
in the x,y direction. These two oscillators share energy due to Larmor precession around the magnetic field (in the 
direction of z" axis). Also, as can be seen from the dipole component of radiation tensor Jq of unpolarized radiation, 
J5 2 oc J7q 2 oc (1 — 1.5 sin 2 r )/2. Therefore the dipole component of the density matrix p\ , which is proportional to Jq 
(Eq[TTJ Wt\ , changes from parallel to perpendicular to the magnetic field at Van Vleck angle (Figj5]) . As the result the 
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Fig. 6. — Degree of polarization of S II absorption lines vs. r , the angle between magnetic field and direction of pumping source (left); 
8, the angle between magnetic field and line of sight (right). 

polarization of the absorption line also changes according to Eq.([T]). This turnoff at the Van Vleck angle is a general 
feature regardless of the specific atomic species as long as the background source is unpolarized and it is in the atomic 
realignment regime. 

The absorption shown in Fig |6l8l ?eft exhibits variations with 9 r and 9. The maximal polarization (around 12 percent) 
is expected for the S II (J u = 1/2) component. At 6 = 90°, the observed polarization reaches a maximum for the 
same 9 r and alignment, which is also expected from Eq.([T]). This shows that atoms are indeed realigned with respect 
to magnetic field so that the intensity difference is maximized parallel and perpendicular to the field (Fig|2]ri<?/i/) . At 
9 = 0°, 180°, the absorption polarization is zero according to Eq.{T]). Physically this is because the precession around 
the magnetic field makes no difference in the x, y direction when the magnetic field is along the line of sight (Fig^fright). 

The optical depth also varies as a result of the alignment (see Eqj3]). We plot the optical depth ratios of line 
(J u = 1/2) to line (J u = 3/2) (see Fig|7]and right panel of Fig[S]). For optically thin lines, the equivalent width ratio 
is equal to that of optical depths. Note that when the magnetic field is parallel to the incident radiation, the magnetic 
field does not have an impact. The corresponding line ratios without magnetic realignment are equal to those values 
at 9 r — 0° in the graphs. 

The fact that absorption polarization can be only parallel or perpendicular to the magnetic field in the plane of sky 
has an important implication. It means that the direction of polarization traces the projection of magnetic field in the 
plane of sky (the angle 4>b, see FigfJJ within a 90° uncertainty. This is a very "easy" information. Furthermore if we 
get two measurables, the degree of polarization and the ratio of a doublet 8 , we shall be able to extrapolate both 9 r 
and 9 = 9b (see FigHJ). With 9 r known, we can remove the 90° degeneracy of 4>b and thus obtain the 3D information 
of the magnetic field (9b,4>b)- 

We consider a general case where the pumping source does not coincide with the object whose absorption is measured. 
If the radiation that we measure is also the radiation that aligns the atoms, the direction of pumping source coincides 
with line of sight, i.e., 9 — 9 r (Figg]). 

Cr I has a ground state 7S3 and excited state 752,3,4- The computing procedure is similar to that for S II. The 
difference is that the density matrix for the Cr I ground state has seven components with k ~ 3 — 3, 3 — 2, ...,3 + 3 = 
0, 1, 2, 3, 4, 5, 6, 7. We only need to consider even components /9q' 2,4 ' 6 . Following the same procedure as described for 
S II, we get the following results: 



V 2 (J*)" 






■( 







35.829 cos 29 r - 12.593 cos 46V + 0.253 cos 66V + 1.932 
.671 cos 26V + 1.2648 cos 40 r - 0.6092 cos 60 r + 1.7263 
-1.0902 cos 261 - 0.5032cos46V - 0.2516 cos 66V - 0.5404, 



/(-322cos26V + 1.24 cos 40 r + 568).(19) 



By inserting these results into Eq.(TT]) or one can get the polarizations of Cr I absorption lines. 

4.2. More complex structure with multiple lower levels: C II, Si II, O I, S I, & Cr II 

4.2.1. Strong Pumping 

In this section we discuss atoms with more complex structure. As illustrated in Table 1, C II has not only mul 
tiple upper levels, but also two levels 2P°, 2 3 ^ 2 on the ground state 



The level 2P°y 2 has altogether two magnetic 



8 S II is a triplet, for which we can get three degrees of polarizations and two line ratios. However, the degrees of polarizations are 
proportional to each other as they have the same dependence on CTq(J;) (see EqfTJl. So do the line ratios. Thus only by combining one 
degree of polarization and one line ratio, we can obtain the two angles. 



11 



Line ratio (J =1/2)/(J =3/2) of S II absorption 



0.6 
0.58 








- - -9=0° 
-9=45° 


0.6 
0.58 


0.56 








- - 9=90° 


0.56 


0.54 


s 

\ 

\ 






/ 

/ 

/ 


0.54 


■C 0.52 


\ 






/ 

/ 


■*£ 0.52 






\ 


~' ~ / 






0.5 




^\ ' 
z 






0.5 


0.48 










0.48 


0.46 










0.46 


0.44 










0.44 



Line ratio (J =1/2)/(J =3/2) of S II absorption 





20 40 60 80 100 120 140 160 180 



20 40 60 80 100 120 140 160 180 



Fig. 7. — Optical depth ratios of S II absorption lines vs. 9 r , the angle between magnetic field and direction of pumping source (left); 
the angle between magnetic field and line of sight (right). 
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Fig. 8. — The contour graphs of S II polarization and line ratio. They are determined by the dipole component of density matrix a"^(8 r ) 
and the direction of observation 9 fEqs lll3] l. Since the 9 r (Fig[T} dependence is solely from Oq, they exhibit similar pattern in the vertical 
direction. 

projections Mj = ±1/2, and is thus not alignable and has only one density tensor p®. Level 2P£/ 2 has density com- 
ponents p"' 2 . We shall consider the five lines 1036.34A(2P° /2 -> 2S 1/2 ), 1037.02l(2P 3 ° /2 -> 2S 1/2 ), 1334.53A(2P° /2 -> 
2D 3/2 ), 1335.66A(2P° /2 -> 2L> 3/2 ), 1335.71 (2P 3 ° /2 -> 2D 5/2 ), which are observed in QSOs. All the transitions between 
upper levels and the ground state should be counted. There is also a magnetic dipole transition between the two ground 
levels. Although its transition probability is very low, it can be comparable to the optical pumping rate in regions far 
from any radiation source. Depending on how far away the radiation source is, there can be two regimes divided by 
the boundary where the magnetic dipole radiation rate A m is equal to the pumping rate t^ 1 . Inside the boundary, 
the optical pumping rate is much larger than the Ml transition rate A m so that we can ignore the magnetic dipole 
radiation as a first order approximation. Further out, the magnetic dipole transition is faster than optical pumping 
so that it can be assumed that most atoms reside in the lowest energy level of the ground state. For the case of C II 
that we consider now, the lowest level is not alignable and so we shall concentrate on the first regime, and neglect the 
magnetic dipole transition. 
The two regimes are demarcated at n (see Fig|9]!e/Z), where r\ is defined by 



t^ 1 =Bj lJu I B B{R*/ri) 2 = A Tl 



(20) 



and Ibb, R* are the intensity and radius of the pumping source. Outside the boundary where r^ 1 <C A m , C II and Si 
II are not alignable and so their absorption lines are unpolarized. For different radiation sources, the distance to the 
boundary differs. For an O type star, the distance would be ^O.lpc, while for a shell star (T e // = 15, 000K), it is as 
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Fig. 9. — left: a carton illustrating how the atomic pumping changes with distance around a radiation source. For circumsteller region, 
magnetic field is stro ng, such that the Hanle effect, which requires ~ A, dominates. Atomic alignment applies to the much more 
distant (< T2, see Eql36|l inte rstellar medium. Inside T2, it can be further divided into two regimes: strong pumping and weak pumping, 
demarcated at ri (see Eq |20l l ; Right: whether and how atoms are aligned depends on their intrinsic properties (transitional probabilities 
and structures) and the physical conditions: density n, temperature T and the averaged radiation intensity from the source 7*. If the 
pumping rate is less than collisional rates, alignment is partially destroyed. Then for atoms with multiple lower levels, depending on 
the comparison between the pumping rate and the magnetic dipole radiation rate among the lower levels, the atoms are aligned differently. 
In the strong pumping case, all th e align able lower levels are aligned; on the contrary, only the ground level can be aligned in the case of 

absence of alignment outside the radius r\ ensures that observations constrain the 
magnetic field topology within this radius. 

Since there are multiple transitions with frequencies that span over a wide range, the optical pumping and its 
resulting atomic alignment depend on the shape of the spectrum of the radiation source. Neglecting the difference 
within an LS multiplct, we characterize the intensity of the multiplet (2P° — ► 2S* 1 / 2 ) by Is, multiplet (2P° — * 2P1/2.3/2) 
by Ip, and that of multiplet (2P° — > 2D 1 / 2 ,3/2,5/2) by Id- In the case that other multiplets of the same terms have 
comparable intensities, Is, Ip, Id are weighted summations of all the intensities of the same terms (e.g, see Eq |24|) . 

Is = A036.3 

Ip = ^904.1 

Id= ^1334.5 (21) 

These are intensities per unit frequency. So are other intensities in the paper. We use wavelengths (in unit of A) rather 
than frequencies in the subscripts as the lines are in optical and UV band. Note that, due to absorption from the 
Lyman series, those lines with wavelengths shorter than 912A can be ignored in the case of distant radiation pumping 
sources. We obtain for C II: 

p°(J, = 3/2) = g° [(o.447^ + 0.39I P Io - 1.65I s Id - 7.067 2 - 7.377 P 7 S ) cos 4 9 r + (-6.547)?, - 47.817 P 7 D - 17.567 S 7 D + 25.8S7 2 
+ 12.797j.7s) cos 2 6 r + 31.337 2 , + 408.707p + 131.637? + 302.237 D 7p + 137.297x57s + 505.487p7 s ] 

p 2 (J, = 3/2) = q° [-0.727 2 - 3.777j=7 D + 8.887 S 7 D + 12.757 2 + 32.917? + (2.157?, + 11.327p7 D - 26.657S7D - 38.257 2 - 98.727? 

- 275.407 P 7 s )cos 2 r + 91.807p7 s ] , (22) 



where 



g° =p°,(Ji = l/2)((-0.607?, + 4.927 P 7 D + 16.207 S 7 D - 9.987 2 - 52.297? - 60.487p7 s ) cos 4 B T + (-4.07?, - 21.137p7 D 

- 23.987 s 7£, + 14.977?, + 34.867? + 45.887p7 s ) cos 2 6 r + 21.957?, + 159. 687?, + 87.157? + 135.087i,7p + 98.797£,7 S + 239.827p7 s ). 



Si II has the same structure as C II does. Though there are many lines in the range from 989. 87A to 1533. 43A, 
their upper levels belong to only three terms: 2Si/ 2 , 2Pi/2,3/2 an d ^Di/2,3/2,5/2, the same as C II. The difference of its 
alignment from C II is solely due to the difference in the line strengths of the multilcvcls involved, 

Po(,/, = 3/2) = Qq [(-0.687?, - 0.167P7U + 0.467 S 7 D + 0.747? + 0.547p7 s ) cos 4 B T + (l0.097?, + 19.327 P 7£, + 4.927 S 7 D - 2.727? 
- 0.937p7 s ) cos 2 6 r - 48.317?, - 42.977? - 6.657? - 121.927£,7 P - 38.367£,7 S - 36.827 P 7 S 1 



p 2 (,/, =3/2) = qI [(7?, + 2. 01937 p 7 d + 0.21967 s 7£, - 0.14127? - 0.17507? - 0.70377 p 7 s ) (9. 497 - 28.490 cos 2 9 r \\ 



(23) 
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3.3543 
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237.11 


342.79 


255.94 


172.05 


Ip/W 
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Id/W 


81.5464 





230.45 





222.49 


76.4647 


218.89 


251.90 


390.29 


314.73 


88.5257 



TABLE 5 

The line intensities (10~ 4 ergs/cm 2 ■ s ■ Hz ■ sr) for pumpings of C II, Si II, O I, S I and Ti II (see EQ I21I24I32I33H in typical 

ASTROPHYSICAL ENVIRONMENT: INTERSTELLAR RADIATION FIELD (ISRF), SHELL STAR (T e ff = 15,000K), O-TYPE STAR (T e ff = 50000K). 

W IS THE DILUTION FACTOR (R»/r) 2 . 

where 

g° = p°(Ji = 1/2)/ [(ll.59/o - 11-96/p/c - 11.73/s-fo + 1.05/J, + 2.63/ 2 + i.40I P I s ") cos 4 0,. + (-0.920/2, + 15. 182 J P J D 
+ 11.5Q6/ S / D - 1.573/p - 1.756/| - 3.338/ P / s ) cos 2 8 r - 32.783/o - 16.781/ 2 , - 4.389/?, 

and 



55.5989/r>/p - 28.387/ D / s - 17.447/p/ s J , 



Is = (-2l533.4<Sl533.4 

Ip = ■'1194.55 

Id = (^1816.9'S'l816.9 



-fl309.3<5i309.3 + -fl023.7Sl023.7)/Sl533.4, 
-^1264. 7^1264. 7 + ■■■)/^'l816.9> 



(24) 



in which S\ oc gB is the line strength of the line with wavelength A. The typical values of these intensities Is, Ip, Id are 
listed in Table If we use these density tensors with Eqs (|ll3[) . we get the absorption coefficients and the corresponding 
polarizations and line ratios. O I and S I have ground states 3Po,i,2 and upper levels 3S1, iD° 2 3 . The ground level 

is 3i-2- Thus, unlike C II and Si II, O I and S I are alignable in both regimes regardless of 1 > or < A m . The degree 
of alignment, however, differs in the two regimes. The case for weak pumping will be presented in ^4.2.21 In the case 
of strong pumping, that we consider here, both levels of 3Pi2 are alignable. Eq. pi[) is then a five dimensional matrix 
with p°' 2 (Ji — 1) and p 0,2,4 (J; = 2). The explicit generic solution is very lengthy. Here as an example, we give the 
results for pumping by a shell star (a black body source of T = 15000K, see Table [5]). For O I: 

p° (Ji = i) 
pliJi = i) 

P°o(Ji =2) 



-P°o 



/0.0851 cos 14 I 
— 1 . 1006 cos 1 
0.1051 cos 14 ( 

Vl.3136 cos 12 ( 



- 0.4823 cos 1 ^ 8 T 
+ 3.8022 cos 10 I 

- 0.1746 cos 12 r 

- 5.3551 cos 10 6 r 



- 1.3183 cos 10 9 r + 7.2195 cos 8 6 
+ 3.4199 cos 8 6 r - 12.6453 cos £ 
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where eg = p"(Ji = 0)/(2.3cos 14 B r - 4.4cos 12 8 r - 57.5cos 10 8 r - 64.3cos 8 8 r + 657.9cos 6 B r + 1620. 9cos 4 8 r - 2519. 7cos 2 8 r - 6083.6). In 

fact, C I, Si I and S III have the same ground state term 3Po,i,2, but with an opposite energy sequence, namely, 3Pq 
is the ground level. Thus they are only alignable in the strong pumping case. In this case, their alignment is similar 
to that of O I and S I. The results are given in App. [F] 

Cr II has a ground state 066*5/2 and excited state z6P^ 2 5 y 2 7 y 2 . Transition between them generates the triplet 2056, 

2062, 2066A. The most important difference between CrII and the other species discussed in this paper is the presence 
of a metastable state a6-Di/2,3/2,...,9/2 (with lifetime tt ~ 10s), which can act as the proxy of the ground state. The 
transitions between the metastable levels and the upper levels are permitted (with wavelength ranging in 2740-2767A). 
An analogy thus can be made between Cr II and the atomic species with multiple lower levels we deal with in this 
paper. In the strong pumping regime, both the ground level and the metastable level can be aligned. Using the similar 
treatment (neglecting the forbidden transition between the metastable levels and the ground level), we can get the 
alignment and polarizations for the two groups of transitions (the triplet from the ground state and the absorptions 
from the metastable state) . Since there are many levels involved and the dimension of the equations we need to solve 
is 17 x 17, we only present selected numerical results here. The dipole components of the density matrices for the 
ground and metastable level is given in F\g®~ight. We show also the polarization of the line (a6S — z6P^ 2 , 2066A) 

and the line (a6D 3/2 - z6P£ /2 , 2750 A) in FigHD] 

There is an advantage to use the absorption lines from the metastable level. Because the energy of the metastable 
level is high, they can be populated only through optical pumping. The absorption lines from the metastable levels 
therefore only traces the medium in the vicinity of a pumping source. Absorption from the ground state, however, 
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Fig. 10. — The contour graphs of Cr II polarization and line ratio. The upper panels are for the absorption from the ground state a6S 
and the lower panels are for the absorption from the metastable state a6D. 

can happen in regions much further away where alignment is marginal. The actual observed signals shall depend on 
N a /Ntoti the ratio of column density to total column density. Combining the absorptions from the metastable state 
and the ground state, allows polarimetric testing of the conditions in the vicinity of the pumping source. Because of 
that we intend to provide a study of more atoms with metastable states (see Table [2]) in our next paper. 

4.2.2. Weak Pumping 

Now let us consider the second regime in which the radiation source is very distant so that the magnetic dipolc 
transition rate within the ground state dominates over the optical pumping rate. Because of the faster magnetic 
dipole decay, we can neglect the optical pumping from levels other than the ground level, i.e., only the ground level is 
alignable. In this case, additional magnetic dipole transitional terms should be added to Eq.©. For the ground level, 
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For those lower levels other than the ground level, 

nk(T,\--X^A ^7\_7'l„,.. \nft( TA _L W_1 \J[ + Jl+k+l\ T >] J J l J l k 



(-l) JMk+1 [J u ] \ J T lJ T lk A A(J U - ,h)p^J u ) (28) 



As discussed earlier, because of fast magnetic mixing, there are no interference terms and thus there are only p§ terms 
in the above equations. Combining with Eq.Q, we get the following set of linear equations: 
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By solving these equations, we can get the ground level (3P2) dipole density component of O I and S I: 

al(J° = 2) = [(5.14D23 2 + 2.21P22D23 + 13.21S21P/23 - 44.74S21 2 + 9.96P22S2l) cos 4 8 r + (-46.47D23 2 - 63.20P22.D23 + 31.17S21D23 (31) 

- 18.59P22 2 + 113.35521 2 + 58.30P22S2l) cos 2 8 r + 14.92D23 2 + 6.20P22 2 - 32.81S21 2 + 20.82D23P22 - 11.86Z323S21 - 20.54P22S2lj / 
[(-2.35P/23 2 + 5.63P22D23 - 5.66S21D23 + 5.45P22 2 + 11.86S21 2 + 2.77P22S2l) cos 4 G r + (21.83P/23 2 + 8.01P22P/23 + 47.35S21D23 

- 22.10P22 2 + 37.97S21 2 + 40.71P22S2l) cos 2 8 r - 71.97P/23 2 - 141.15P22 2 - 133.61S21 2 - 215.84D23P22 - 199.53.D23S21 - 281.43P2252lj 

where S21, P22, D23 arc = £ I V S V . For O I, 

5*21 = Is5i302.2 = 1-08/5 = ^1302. 25*1302. 2 +^1039.2 5*1039.2, -P22 = 0, 

D23 = I]jSgs8.7 = 0.732/d = /1025 5*1025 + ^988. 7^988. 7 + ^971. 7^971. 7. (32) 

for S I, 

521 = IsS*i826.2 = 0.502ii 8 26.2, -P22 = ip5*i295.6 = l-78/i295.6j 

D23 = /d5i425 = 0.149/d = Zl4745i474 + ^14255i425 + ^1316. 55*1316.5 (33) 

The line strengths S v in the above two equations are taken from NIST Atomic Spectra Database. The line intensities 
for various environments are listed in Table 5, where the values for interstellar radiation field (ISRF)are particularly 
relevent to the weak pumping regime we discuss here. 

We only provide the density matrices for the ground level. Absorption from the levels other than the ground level 
is negligible because the absorption rate is less then the magnetic decay rate. 

4.2.3. Role of collisions 

In the above calculations, collisions are neglected. Now let us see whether this is a reasonable assumption. Basically 
there are two types of collisions. The first type is an inelastic collision, which excites transitions among the J levels of 
the ground state. The second type of collision is elastic and only causes the angular momentum J to flip. 

The inelastic collision can happen with electrons or hydrogen atoms. Both rates can be found in Spitzer (1978). The 
transition rate for collisions with electrons is 

8.63 x 10- 6 n e q(j, fc) 3 _ x 
fhj = n e <y k j = cm 8 ' ( 34 ) 
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where [k] is the statistical weight of level k. a, the collisional strength, is of order unity (Osterbrock 1974), while 
the magnetic dipole emission rate between the lower levels is A m > 10 -6 s -1 . Apparently, the collisional rate is much 
smaller for diffuse medium. For collisions with hydrogen atoms, fkj — nn^kj ^ 10 _10 7ii7S _1 ; this is also smaller than 
the magnetic dipole decay rate A m . Thus both collisions with electrons and hydrogen atoms are less frequent than 
magnetic dipole decay in diffuse medium. 

The second collision involves spin flip and is dominated by Van der Waals interactions with neutral atoms, in 
particular, hydrogen atoms. The corresponding rate is (Lam & ter Haar 1971, Landi Degl'Innocenti 2004), 



fsf - 2 x 10 



-10 



13.6eV 

Ei — Ei 



0.8 



T(l + ±) 



0.3 



n H s \ (35) 



where Ei is the ionization potential, Ej is the energy of J level, p is the atomic weight. 

In the weak pumping regime which we discussed in £ )4.2.2[ collisions can become comparable to the optical pumping 
rate. In this case alignment is partially destroyed 9 . Since pumping rates for different species are different, we can 
find the tomography of the magnetic field by utilizing lines from various species. By equating the collisional rate 
max(/fcj, f s f) to the optical pumping rate, we obtain the alignable range for an O type star, 

for C II 

(36) 




For a typical QSO with luminosity < L >= XL\ = 10 46 16 ergs • s^ 1 at HOOAand a spectrum (Telfer et al. 2002) 
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the alignable range would be 



1L H nT VuFkJ Mpc,forSiII 

In actual observations, there is inevitably an average along line of sight. The absorption from atoms outside the 
alignable range tend to reduce the degree of polarization. This adds another dependence on N a /N tot , the ratio of 
alignable column density to total column density. In this case, there is an advantage to choose species present only 
near a pumping source. For instance, the column densities of many ions within an HII region may dominate the 
column densities along the line of sight. In this case, the conditions in the HII region will be traced. 

5. SUMMARY OF THE PROCEDURE ADOPTED 

In this paper, we provided calculations for the alignment of atoms with fine structure owing to pumping by anisotropic 
radiation in a magnetized medium. The angular momentum of the radiation is transferred to the atoms and causes 
differential occupation in terms of angular momentum (magnetic quantum number) of the ground state. This results 
in atomic alignment. Any background radiation (including the pumping source) coming through these atoms will be 
absorbed differentially, causing polarization in absorption lines. 

In most astrophysical environments, incident radiation is axially symmetric. This makes the irreducible density 
matrix the suitable formalism to describe the system (see App(XJ. In this case, alignment is simply represented by 
the p 2 q term, which causes linear polarization. Total population is given by the p% term. Since the incident radiation 
is unpolarized, there is no orientation, and the corresponding terms p x q are zeros accordingly. Circular polarization, 
resulting from orientation, does not exist in this case 10 . 

When the magnetic field exists, it is convenient to choose the magnetic field direction as the quantization axis. 
Coherence terms Pq^q — disappear because of fast magnetic mixing on the ground state. In this sense, atoms are 
realigned along the magnetic field compared to the case without a magnetic field. This requires a transformation from 
the general observational frame to the theoretical frame where the magnetic field defines the z axis (Fig[3|. This can 
be done by two Euler rotations (see App |E|) . 

We solve the statistical equilibrium equations of both ground state and excited state and obtain the density matrix 
Pq of the ground states. Multilevel atoms are considered as well. All of the transitions between the levels in ground 
state and the levels in upper state should be taken into account to get the correct atomic alignment. If there are 
multiple lower levels and the decay rate between them is larger than the optical pumping rate, the transitions between 
them should be also taken into account and their corresponding statistical equilibrium equations should be included 

9 To calculate the residual alignment, one should include collisional transitions in the statistical equilibrium equations of the atomic 
populations. 

10 If the background source is polarized in a different direction with alignment, circular polarization can be generated due to dephasing. 
This is a second order effect, and we shall discuss it in our companion paper. 
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to solve the level populations. These transitions can substantially change the alignment. This is the weak pumping 
case we studied in §4.2.21 

Then we insert the dipole component of the density matrix of the ground state gg in to the equation for absorption 
coefficients fEq fMl) . yielding line ratios and polarizations for absorption lines (Eq dll3p . Polarization of absorption is 
either parallel or perpendicular to magnetic field in the pictorial plane because of fast precession around the magnetic 
field in the ground state. 

Fig^right illustrates when and in what form we expect to see atomic alignment depending on the atoms and their 
environmental conditions. It also illustrates the computational steps to determine the direction of magnetic field. 

6. DISCUSSION 
6.1. Choice of atoms 

We have chosen a particular set of atoms (see Table [T]) as examples of alignment for atoms with fine structures. 
Since the only requirement for atomic alignment is having a ground level with total angular momentum J > 1/2, there 
are many more atomic species that can be explored. Regardless of the differences in specific atomic structures, the 
physics is similar. By considering transitions from multiple levels to multiple levels, we have developed a procedure 
that can be applied to virtually all atoms fine structure. We did not include level coherence as it is negligible except for 
a highly turbulent medium where Sv > lOOkm/s. For multiplets, the absorptions Ji > J/ are usually more polarized 
than other components which provides more routes and more balance between a and 7r transitions. 

Practically, however, there is inevitably an averaging along line of sight, which adds another dependence on the 
ratio of alignable column density and total column density N a /N tot . For different species, this ratio is different. The 
ratio should be close to 1 for highly ionized species which only exist near radiation sources. The same is true for the 
absorption from metastable state, for instance, for CrII that we discussed earlier. The energies between the upper and 
metastable state in this case are too high to be affected by collisions with thermal particles in cool phases along the 
line of sight. Thus the corresponding absorption lines are also mostly affected by the medium in the vicinity of a star 
(or QSOs). Combining different species (with different N a /N to t), it is possible to get a tomography of the magnetic 
field in situ. To extend the technique we shall present elsewhere calculations for more atoms with metastable states. 

The particular choice of atoms to use depends both on the instruments available and the object to be studied. 
With the advent of space-based UV polarimetry many more alignable transitions will become important. This will 
allow detailed studies of magnetic fields in circumstellar regions, diffuse ISM, intergalactic medium, supernovae rem- 
nants 11 , QSOs etc. However, even at present a substantial progress may be achieved using ground-based and possibly 
sounding rocket observations. 

6.2. Wavelengths employed and particular cases 

The incomplete list of astrophysically interesting transitions in Table 1 1 121 shows that many of the transitions are in 
the UV part of the spectrum. At the moment, the optical polarization is easily available, e.g., Cr I, Ti II. Furthermore, 
we expect that in the years to come, UV polarimetry will become an essential tool for astrophysical research. In the 
mean time, a number of transitions (Ti II, Cr I, Fe I, etc.) as well as atomic transitions shifted due to cosmological 
rcdshifts can be studied using present ground based polarimetry facilities. The latter is applicable to, e.g., Cr II 
detectable in QSO absorption (Kulkarni et al. 2005). 

We did not discuss transitions beyond the Lyman limit 912A. There are actually many lines in this range, especially 
for absorptions in QSOs (Verner, Barthel & Tytler 1994) and Lyman a clouds (see Akerman et al. 2005). Indeed for 
diffuse ISM, absorptions are severe. However, for QSOs, this constraint is relieved because of proximity effect, which 
is the measured decrease in the number of Lya clouds induced by the UV radiation from a QSO in its vicinity (Madau 
1992). When the light reaches those Lya clouds far apart, it is already red-shifted to a wavelength longer than 912A 
so that it is free of the absorptions from H I as well. Similarly these lines are also informative for the local bubble (see 
Redfield & Linsky 2004). 

As we discussed for the case of weak pumping f i)4.2.2[) . the forbidden decay of the ground state may result in the 
loss of alignment for atoms like C II and Si II, which is an interesting observable effect in itself. For other transitions 
(e.g., O I, S I), the decay does not erase but modifies the alignment. 

It is worth mentioning that when the photon arrival rate is comparable the Larmor frequency, the polarization 
becomes sensitive to the amplitude of magnetic field. This opens an avenue for obtaining the strength of weak magnetic 
fields, especially for circumstellar region and QSOs. If there are several species spatially correlated, it is possible to 
measure both the magnetic field strength and direction. In addition, even measurements of unpolarized lines are useful 
in conjunction with studies of aligned species. If the photon arrival rate is less than Larmor precession rate, atoms are 
aligned with respect to the radiation field and their line polarization does not contain information about the magnetic 
field. For instance as we pointed out in £12. 2[ magnetic field can only realign CII at a distance > 7.7Au from an O star 
if the magnetic field strength ~ 5/iG. 

We would like to emphasize here that we deal with optically thin case, or unsaturated absorption in this paper, 
radiative transfer needs to be considered for the saturated regime. 

6.3. Utility of transitions 
11 This possibility was mentioned to us by Ken Nordsicck. 
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There are a number of parameters that determine the polarization: the direction of the magnetic field and radiation 
field (0, 9 r ), the percentage of radiation anisotropy W a /W and the ratio of alignable column density and total column 
density N a /N to t- The list of observables we have are the degree of polarization and the line intensity ratios (equal to 
optical depth ratios for optically thin case) for multiple lines. Since there are a lot of alignable atomic species, we can 
have more known observables to determine the magnetic field if we make sure that these atomic species coexist. One 
can also use the transitions between different levels of the same atom. 

Combining different atomic species, we can make a tomography of magnetic field as the ratio N a /N tot varies with each 
atomic species. Moreover, the fact that the radius of alignment changes for different environments, allows polarimetric 
testing of the conditions in the vicinity of the pumping source. 

If the alignment is produced by the background radiation source which is used to study the absorption lines 12 (see 
Fig. 1, right), the line of sight coincides with the direction of radiation, resulting in r = 0. 

Our calculations in this paper were aimed at exploring the feasibility of the measurements of the astrophysical 
magnetic fields using atomic alignment. In the paper above we dealt with the fine structure line transitions, while the 
hyperfine structure transitions are considered in a separate paper. Taken together, the alignable species allow studies 
of magnetic fields and environmental conditions for various astrophysical systems, from comet atmospheres to stellar 
winds, supernovae, reflection nebulae, interstellar gas and the intergalactic medium. 

6.4. Studies using contour plots 

Practically, we can make use of contour plots like what we did for S II, Cr II to get information on the magnetic field 
and environment from observational data. For instance, if one measures for S II (J u — 1/2) absorption line degree of 
polarization per unit optical depth 8% and the equivalent width ratio of S II (J u = 1/2,3/2) doublet 0.48 13 . First 
from the direction of polarization, we know that the magnetic field projection is either parallel or perpendicular to the 
direction of polarization in pictorial plane (azimuthal angle ips, see Fig|4]). This is already the whole information one 
can expect from the Goldreich-Kylafis effect (Goldreich & Kylafis 1982). Then, using the degree of polarization and the 
line ratio, one can get further information about the magnetic field in the plane of line of sight 0, and also the direction 
of pumping source r - By overlapping the two contour graphs (Fig[5]), we can find their crosses at r ~ 90° ± 72° and 
~ 90° ± 7° 14 . With Or known, we can easily remove the 90° degree degeneracy of 4>b according to the Van-Vleck 
effect and thus obtain the 3D information of the magnetic field. If we know two species coexist spatially, we can also 
use their ratio of degrees of polarizations and their ratio of equivalent width. In this case, we don't need to know the 
optical depth for individual line. Alternatively, r and can be obtained by solving Eqs (|ll3|) . Calculating r and 
several ways cross-check the results and decreases the uncertainties. 

6.5. Complementary techniques 

Atomic alignment can render unique information about the magnetic field. For instance, grain alignment 15 (see 
review by Lazarian 2003) provides information about the component in the plane of sky. Similar to the case of 
interstellar dust, the rapid precession of atoms in a magnetic field makes the direction of polarization sensitive to 
the direction of the underlying magnetic field. As the precession of magnetic moments of atoms is much faster than 
the precession of magnetic moments of grains, atoms can reflect much more rapid variations of magnetic field. The 
limited radius of atomic alignment is also an advantage to enable local measurement of 3D magnetic fields while grain 
alignment tests global 2D topology. More importantly, alignable atoms and ions can reflect magnetic fields in the 
environments where either the properties of dust change or the dust cannot survive. This opens wide avenues for 
magnetic field research in circumstellar regions, interstellar medium, interplanetary medium, intracluster medium and 
quasars, etc. 

The complementary nature of atomic alignment to the Goldreich-Kylafis effect 16 (Goldreich & Kylafis 1982; Girart, 
Crutcher, & Rao 1999). is that the former is applicable to the regions near emitting sources (e.g. circumstellar regions), 
while the latter is applicable to molecular clouds. Moreover, atomic alignment provides 3D tomography of magnetic 
fields that is not measurable otherwise. 

A detailed discussion of alignment of atoms in different conditions corresponding to circumstellar regions (see Diner- 
stein, Sterling & Bowers 2006), AGNs (see Kriss 2006), Lyman alpha clouds (see Akerman et al. 2005), interplanetary 
space (see Cremonese et al. 2002), the local bubble (see Redfield & Linsky 2002, 2004), etc., is given in our companion 
paper. Note, that for interplanetary studies, one can investigate not only spatial, but also temporal variations of 
magnetic fields. This can allow cost effective way of studying interplanetary magnetic turbulence at different scales. 

The change of the optical depth is another important consequence of atomic alignment. The actual calculations 
should take into account the realignment caused by a magnetic field. It may happen that the variations of the optical 

12 We call this case a "degenerate" case to distinguish from a case where the background source and the source of the aligning radiation 
differ. 

13 For optically thin lines, the equivalent width ratio is equal to optical depth ratio as given in the lower panels of Fig|8] 

14 The uncertainties come from the degeneracies of sin 2 9 and cos 2 9 r in Eqs llll3jl . 

15 Grains have tendency to be aligned with their long axes perpendicular to the magnetic field. Curiously enough, the interaction with 
photons is also important for grain alignment. 

16 Atomic alignment has some similarity to Goldreich-Kylafis effect, which also measures magnetic field through magnetic mixing. 
However, Goldreich-Kylafis effect deals with the polarization of radio lines. The upper states of radio lines are so long lived that significant 
magnetic mixing can happen among different magnetic sublevels of these states. Atomic alignment, on the other hand, happens with 
ground states of optical and UV transitions. 
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depth caused by alignment can be related to the Tiny-Scale Atomic Structures (TSAS) observed in different phases of 
interstellar gas (see Heiles 1997). 

7. SUMMARY 

In this paper we identified a set of atoms and ions with fine structure that provide substantially intensive absorption 
lines in the visible and near UV wavelength range. We calculated the degree of alignment of these species and obtained 
the relation between the Stokes parameters of the absorbed radiation and the 3D orientation of magnetic field with 
respect to the observer and the pumping/absorbed sources. Our calculations can easily be extended to other absorbing 
species with fine structure. In particular, we have shown that: 

1. The ground-state alignment happens as the result of interaction of atoms and ions that have > 2 sublevels on the 
ground state with the anisotropic flow of photons. If the rate of excitation from the ground state is lower than the 
rate of atomic precession in the external magnetic field, the alignment happens with respect to the magnetic field. 

2. The ground-state alignment affects the polarization state of the absorbed light as well as the intensity of different 
line components. Combining information from different components and absorption lines one can study both 3D 
topology of magnetic fields and physical conditions in important astrophysical environments that include quasars, 
interstellar medium and extragalactic gas, etc. 

3. For the light absorbed by aligned atoms, the direction of polarization is either parallel or perpendicular to the 
magnetic field in the plane of sky. The switch between the two possibilities happens at the Van-Vleck angle between 
the direction of the magnetic field and the pumping radiation. 
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APPENDIX 
IRREDUCIBLE DENSITY MATRIX 

In quantum mechanics, a particle (atom or photon) can be described by its wave function \ip >= ^2j M cjm\JM >. 
An equivalent but more elegant description without the phase arbitrariness is a density matrix p — \ip >< ip\ = 
^Zj M ji M i cjmc*j, m ,\JM >< J'M'\. The diagonal term JM\p\JM >= \cjm\ 2 represents the probability of finding 
the particle on state \ JM >; off-diagonal terms < JM\p\ J'M 1 >= cj M c* JIM , characterize the coherence between two 
different states \JM > and | J'M' >. The advantage of the description by density matrix becomes apparent when we 
are dealing with an ensemble of n particles without mutual interactions. In this incident, the density matrix of the 
ensemble is simply the arithmetic mean of individual densities — Y17=i 

When magnetic field is present, it is convenient to choose the magnetic field direction as the quantization axis for an 
atom because Mj remains a good quantum number on the ground state. Meanwhile, the density matrix of radiation 
(or their classical equivalent Stokes parameters, see App[D]) should be obtained in a basis where the line of sight is the 
quantization axis. This requires several rotations of the successive reference frames (see App|E]). Usually, the incident 
radiation is axially symmetric, while the spontaneous emission is spherically symmetric. This makes the irreducible 
tensorial operators pg the suitable formalism for performing the calculations (Fano 1957, Dyakanov & Perel 1965, 
Omont 1965, Bommicr & Sahal-Brechot 1978). The relation between irreducible tensor and the standard density 
matrix is 

Pq W J') = E (-l) J ~ M (2^ + 1) 1/2 ( J M g J M , ) < JM\p\J'M> > . (Al) 

MM' \ / 

For photons, we shall consider only their electric dipole part, which has a spin J = 1. The irreducible tensor has 
a clear physical meaning. In particular, p[j is proportional to the overall population of an atomic level or the total 
intensity of the radiation; p^ describes orientation of an atom or circular polarization of radiation; pg represents the 
alignment term, which are responsible for linear polarization; the terms with Q ^ are coherence terms. For atoms, 
Pq (J, J') is the level-crossing coherence (Bommier & Sahal-Brechot 1978), Pq^ (J) is the Zeeman coherence as they 
correspond to the crossing term pmm> according to the selection rule of 3j symbol (App lCj) in Eq. (|Aip . Note that in 
the regime we study (namely, weak magnetic field [< 1G] and low temperature diffuse medium), individual J levels 
are much narrower than their separations so that level-crossing coherence is negligible. 

RADIATIVE TRANSITIONS 

We consider atoms with fine structure. The ground level of the atoms has an angular momentum J, M and the 
excited state is J', M' 17 . For spontaneous emission from a state J'M 1 of higher energy to a lower energy state J, M, 



' Throughout this paper, we follow the convention that M represents the projection of angular momentum along magnetic field— the 
quantization axis. 
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the transition probability per unit time is 



_64^^ E|<jM|n|j , M , >|2 (m) 

Q 

where V q '° — r • e*'° is the projection of the dipolc moment along the basis vector e q of radiation, e± = (^x- y)/\/2, 
e = z. According to the Wigner-Eckart theorem (Weissbluth 1978), the electrical dipole matrix element is 

Rj ,j =< JM\V q \J'M< > = (-if+M-i (j^',) < J\\V q \\J> >, (B2) 

where the matrix with big " ( )" is Wigner 3j symbol, which is a measurement of coupling of ground angular momentum 
(J, M) and photon angular momentum (l,q) forming angular momentum (J',M') of the upper level (see App. A). 
< J||tq|| J' > is the square root of the line strength, whose value can be found in atomic physics books, e.g., Cowan 
(1981). 

COUPLING OF ANGULAR MOMENTUM VECTORS 

Each atomic (or photon) state has a particular angular momentum associated with it. Because of its vector nature, 
we can describe it with two quantum numbers \jm >, where j is total angular momentum and m is its projection along 
the quantization axis. Whenever discussing any atomic processes involving direction, we would find it is inevitable to 
involve angles and momenta. Angular momentum theory arises as a result. 

When two angular momenta J, mi and J',m 2 are coupled, the resultant angular momentum j,m can be described 
as (Zare 1988) 



\jm >= < Jrai, J'm 2 \jm > | Jm\, J'm 2 >, (CI) 



E 

mi ,m2 



where the elements of transformation < Jmi, J'm 2 \jm > are called Clebsch-Gorden coefficients or Wigner coefficients. 
To better reveal its symmetric properties, the Wigner 3-j symbol was introduced 



J J' .13 

mi m 2 — m 3 



) = (-l) J - J '- m3 b'3]- 1/2 < Jm u J'm 2 \jm > 



(C2) 



The three J must satisfy the triangle relations: 

J + J' > J3, J' + .13 >J,J + 33> J'- (C3) 

and mi + m 2 — m 3 = 0. If one is concerned radiative transitions, one of the angular momentum (J',m 2 ) would be 
that of photon (l,q). It then follows that the selection rules for transitions between states J, mi, jz,mz are 

Aj = j 3 - J = 0, ±1, Am = — m\ = 0, ±1 (C4) 

with the restriction that J = j'3 = is not allowed. 

When three angular momenta J, mi J', mi j'3, need to be coupled, there is more than one way to do this. We can 
first couple J, mi and J', m 2 to get ji 2 , then add j'3; or we can first couple J3, m.3 and J', m 2 and get the same result. 
The transformation between < ji2j3jm > and < Jj 2 3 > is given by (Zare 1988) 

\Jj23jm >= < h2hj\Jj233 > \h2j3jm > (C5) 
312 

where the expansion coefficients < ji2j3j |</j23j > are called rccoupling coefficients as they describe the transformation 
between two coupling schemes < ji2j3jm > and < J) 2 3 >• A closely related quantity was then defined by Wigner 
(1965), which has higher symmetry: 

( / J - f ) = (-l) J+J ' +h+j \ji2,h3}- 1/2 < 3W33\Jj233 > (C6) 
L 33 3 323 J 

Similarly, the recoupling of four angular momenta leads to the 9-j symbol: 



J J' 312 
ji 33 334 
.714 J23 j 



[jl2,j23,j34, J14] 1/2 < jl2j34j|jl4j23j > (C7) 



STOKES PARAMETERS 



The electric field in the plane perpendicular to the propagation direction can be developed over the usual basis 

( e a;j e y)- 

E(t) = [£ie x (cosu;£ - 5i) + E 2 e v {cosujt - 6 2 )] (Dl) 
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When 5i = S 2 , the light is linearly polarized. For theoretical calculations, it is more convenient to use basis composed 
of left- (e_) and right- (e+) handed circular vibrations. They are defined as 



e± = [±ei 



ie 2 



}/V2. 



By using (e + , e_), one has a radiation tensor 



n 



< E+E+ > < E+E- > 

< E-E+ > < E-E^ > 



n++ n+_ 
n_ + n 



(D2) 
(D3) 



where <> denotes an average over time. 
The definition of Stokes parameters Si — [I, Q, U, V] (i=0, 1, 2, 3) and their relation with II is as follows: 

j= J E 1 2 + £ 2 2 = n ++ + n__ ) 
Q = £ 2 -£ 2 2 = -(n + _+n_ + ), 
U = 2E 1 E 2 cos(5i - 82) = -*(n+- - n_+), 

V = 2E 1 E 2 sm(S 1 - 5 2 ) = U++ - 

The observable quantities: intensity, degree of polarization p, direction of linear polarization \ (Fig^right) have simple 
relations with Stokes parameters: 



(D4) 



I = Tr(IL), p=VQ 2 + U 2 /I, sin(2 X ) = U/I, cos(2 X ) = Q/L 
FROM OBSERVATIONAL FRAME TO MAGNETIC FRAME 



(D5) 



In real observations, the line of sight is fixed, and the direction of the magnetic field is unknown. Thus a transforma- 
tion is needed from the observational frame to the theoretical frame where the magnetic field is the quantization axis. 
This can be done by two Euler rotations as illustrated in FigHl In the original observational coordinate (xyz) system, 
the direction of radiation is defined as the z axis, and the direction of magnetic field is characterized by polar angles 9 r 
and 4>b- First we rotate the whole system by an angle <f>B about the z-axis, so as to form the second coordinate system 
x'y'z'. The second rotation is from the z(z') axis to the z" axis by an angle 9 r about the y'-axis. Mathematically, the 
two rotations can be fulfilled by multiplying rotation matrices, 



cos 9 r — sin 9 r 

1 
sin 9 r cos 9 r 



cos 4>b sin (f>B 
— sin <j> b cos 4>b 
1 



cos 9 r cos <fi b cos 9 r sin <pB — sin 9 r 
— sin d>B COS d>B 



sin 



— sin 4>b 

cos 4>b sin 



cos q>B 
~ sin (fig cos 9 r 



(El) 



MORE RESULTS OF CI, SI I, S III FOR CIRCUMSTELLAR CASE 

For a T=15000K black body pumping source, we have obtained the following results, 
C I: 



p?,(Ji = i) 

P 2 o(Ji = 1) 
P°o(Ji =2) 
pliJi =2) 



/ 3.1 cos 10 6 r - 50. 8 cos 8 6 r + 327.1 cos 6 6 r - 
0.01 cos 12 9 r - 0.44 cos 10 6 r - 2.39 cos 8 9 r 
cos 10 6 r - 14 cos 8 e r + 65 cos 6 6 r - 30 cos 4 

\ -0.4 cos 10 6 r - 15.7 cos 8 r + 201.4 cos 6 S, 



1043. 7cos 4 6 r - 1839.1 cos 2 6 r + 7806.8 \ 
-13.13 cos 6 9 r - 29.63 cos 4 6 r - 1367.38 cos 2 9 r + 459.60 
9 r - 3140 cos 2 9 r + 10195 

- 1430.7 cos 4 e r + 4352.5 cos 2 6 r - 1299.1 J 



where = Po(Ji = 0)/(-0.7cos 10 6 r + 11.5 cos 8 6 r - 86.0 cos 6 9 r + 548.9 cos 4 6 r - 1746.6 cos 2 6 r + 4625.5) 
Si I: 



pI(Ji = i) . 

P°(Ji = 2) " 
PliJi =2) 

where qI = p%(J t 
S III: 



So. 



(Fl) 



/ -0.4 cos 10 9 r + 7.9 cos 8 8 r - 62.9 cos 6 8 r + 211.3 cos 4 8 r + 363.1 cos 2 8 r - 1516.6 
-0.013 cos 12 8 r + 0.451 cos 10 8 r - 2.671 cos 8 B r + 18.587 cos 6 B r + 8.482 cos 4 8 r + 150.334 cos 2 8 r - 51.711 
-0.1 cos 10 8 r + 1.9 cos 8 8 T - 4.9 cos 6 B T - 18.1 cos 4 8 r + 635.4 cos 2 B T - 1975.1 

V-0.02 cos 12 8 T + 0.36 cos 10 8 T + 0.38 cos 8 8 r - 23.62 cos 6 8 r + 281.19 cos 4 8 r - 948.70 cos 2 8 r + 285.86 / 

= 0)/(-0.05 cos 12 9 r +0.38 cos 10 6> r -3.23 cos 8 8 r + 10.97 cos 6 r -79.56 cos 4 6> r + 89.78 cos 2 6» r -895. 



(F2) 



20 cos 10 8 r 



1 52 cos 10 8 T + 30 cos 8 8 r - 5710 cos 6 8 r + 23040 cos 4 8 r + 29549 cos 2 8 r - 130181 \ 

- 544 cos 8 B r + 3562 cos 6 8 T + 1663 cos 4 8 r - 393 cos 2 B r - 179 Q„, 
-25 cos 10 8 r + 187 cos 8 8 r - 590 cos 6 8 r - 1259 cos 4 8 r + 54726 cos 2 8 r - 165422 
V 19 cos 10 6 r - 416 cos 8 9 r + 1220 cos 6 6 r + 21954 cos 4 8 r - 91329 cos 2 8 r + 27964 / 

where gg = pg(J ; = 0)/(-3cos 12 9 r + 93 cos 10 9 r - 625cos 8 6» r - 425cos 6 6» r - 155 cos 4 r + 25678 cos 2 9 r - 78043). 



Poi-h = 1) 

pl{Jl = 1) 

P§W =2) 
PoW = 2) 



(F3) 
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